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Abstract. 

The orbital evolution of more than 22000 Jupiter-crossing objects under the 
gravitational influence of planets was investigated. We found that the mean collision 
probabilities of Jupiter-crossing objects (from initial orbits close to the orbit of a 
comet) with the terrestrial planets can differ by more than two orders of magnitude 
for different comets. For initial orbital elements close to those of some comets (e.g. 2P 
and 10P), about 0.1% of objects got Earth-crossing orbits with semi-major axes a<2 
AU and moved in such orbits for more than a Myr (up to tens or even hundreds of 
Myrs). Results of our runs testify in favor of at least one of these conclusions: 1) the 
portion of 1-km former trans-Neptunian objects (TNOs) among near-Earth objects 
(NEOs) can exceed several tens of percents, 2) the number of TNOs migrating inside 
solar system could be smaller by a factor of several than it was earlier considered, 3) 
most of 1-km former TNOs that had got NEO orbits disintegrated into mini-comets 
and dust during a smaller part of their dynamical lifetimes if these lifetimes are not 
small. 
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1. Introduction 

Trans-Neptunian objects (TNOs) are considered to be one of the main 
sources of near-Earth objects (NEOs). Bottke et al. (2002), Binzel et 
al. (2002), and Weissman et al. (2002) believe that asteroids are the 
main source of NEOs. Duncan et al. (1995) and Kuchner et al. (2002) 
investigated the migration of TNOs to Neptune's orbit, and Levison 
and Duncan (1997) studied their migration from Neptune's orbit to 
Jupiter's orbit. Based on the results of migration of Jupiter-crossing 
objects (JCOs) with initial orbits close to the orbit of Comet P/1996 
R2 obtained by Ipatov and Hahn (1999), Ipatov (1999, 2001) found 
that 10-20% or more of the 1-km Earth-crossers could have come from 
the Edgeworth-Kuiper belt into Jupiter-crossing orbits. In the present 
paper we consider a larger number of JCOs than before. Some prelim- 
inary results were presented by Ipatov (2002a-b), who also discussed 
the formation of TNOs and asteroids. The results of the runs of JCOs, 
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including several figures, can be also found in (Ipatov and Mather, 
2003a-b). A wider review on the migration of asteroids and comets to 
NEO orbits was made by Ipatov (2001). 



2. Migration of Jupiter- Family Comets to the Terrestrial 

Planets 

As the migration of TNOs to Jupiter's orbit was investigated by several 
authors, we have made a series of simulations of the orbital evolution 
of JCOs under the gravitational influence of planets. We omitted the 
influence of Mercury (except for Comet 2P/Encke) and Pluto. The 
orbital evolution of more than 9000 and 13000 JCOs with initial pe- 
riods P a <20 yr was integrated with the use of the Bulirsch-Stoer and 
symplectic methods (BULSTO and RMVS3 codes), respectively. We 
used the integration package of Levison and Duncan (1994). 

In the first series of runs (denoted as nl) we calculated the evolution 
of 3100 JCOs moving in initial orbits close to those of 20 real comets 
with period 5<P a <9 yr, and in the second series of runs (denoted as 
n2) we considered 10000 JCOs moving in initial orbits close to those of 
10 real comets (with numbers 77, 81, 82, 88, 90, 94, 96, 97, 110, 113) 
with period 5<P a <15 yr. In other series of runs, initial orbits were 
close to those of a single comet (2P/Encke, 9P/Tempel 1, lOP/Tempel 
2, 22P/Kopff, 28P/Neujmin 1, 39P/Oterma, or 44P/Reinmuth 2). In 
order to compare the orbital evolution of comets and asteroids, we 
also investigated the orbital evolution of asteroids initially moving in 
the 3:1 and 5:2 resonances with Jupiter. For JCOs we varied only the 
initial mean anomaly v in an interval less than several tens of degrees. 
Usually in one run, there were 250 JCOs or 144 asteroids. For asteroids, 
we varied initial values of v and the longitude of the ascending node 
from to 360°. The approximate values of initial semi-major axes, 
eccentricities and inclinations of considered comets (a , e , and i ) are 
presented in Table I. We investigated the orbital evolution during the 
dynamical lifetimes of objects (at least until all the objects reached 
perihelion distance q>6 AU). 

In our runs, planets were considered as material points, so literal 
collisions did not occur. However, based on the orbital elements sampled 
with a 500 yr step, we calculated the mean probability P of collisions. 
We define P as P^/N, where Ps is the total probability of collisions of 
N objects with a planet during their lifetimes, the mean time T=Ts/iV 
during which perihelion distance q of an object was less than the semi- 
major axis a p i of the planet, the mean time spent in orbits with 
aphelion distance Q<4.2 AU, and the mean time Tj during which an 
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Table I. Semi- major axes (in AU), eccentricities and inclinations of several 
considered comets 
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2P/Encke 


2.22 


0.85 


11.7° 


9P/Tempel 1 


3.12 


0.52 


10.5° 


lOP/Tempel 2 


3.10 


0.53 


12.0° 


22P/Kopff 


3.47 


0.54 


4.7° 


28P/Neujmin 1 


6.91 


0.78 


14.2° 


39P/Otcrma 


7.25 


0.25 


1.9° 


44P/Reinmuth 2 


3.53 


0.46 


7.0° 


88P/Howell 


3.13 


0.56 


4.4° 


96P/Machholz 1 


3.04 


0.96 


60.2° 


113P/Spitaler 


3.69 


0.42 


5.8° 



object moved in Jupiter-crossing orbits. The values of P r =10 6 P, Tj, T^, 
and T are shown in Table II. Here r is the ratio of the total time interval 
when orbits are of Apollo type (a>l AU, q< 1.017 AU) at e<0.999 to 
that of Amor type (1.017<g<1.3 AU) and T C =T/P (in Gyr). In almost 
all runs T was equal to the mean time in orbits which cross the orbit 
of the planet and 1/T C was a probability of a collision per year. 

In Table II we present the results obtained by the Bulirsch-Stoer 
method with the integration step error less than ee[lCP 9 -l(U 8 ] and 
also with £<10~ 12 and by a symplectic method with an integration 
step d s <W days. For these three series of runs, the results obtained 
were similar (except for probabilities of close encounters with the Sun 
when they were high). For d s =30 days for most of the objects we found 
similar results, but we found a larger portion of the objects that reached 
Earth-crossing orbits with a<2 AU for several tens of Myr and even 
inner-Earth orbits (IEOs, i.e. with Q<0.983 AU). These few bodies 
increased the mean values of P by a factor of more than 10, and the 
mean probabilities were greater than for d s < 10 days. 

The results can differ considerably depending on the initial orbits of 
comets. The values of P for Earth were about (l-4)xl0~ 6 for Comets 
9P, 22P, 28P, and 39P. For Comet 10P they were greater by an order of 
magnitude than for 9P, though initial orbits of 9P and 10P were close. 
This is a real difference in dynamics of two comets and is not "luck of 
the draw" in the integrations. P exceeded 10~ 4 for Comet 2P. 

The probability of a collision with Earth (or with Venus and Mars) 
for one object that orbited for several Myr with Q<4.2 AU could be 
much greater than the total probability for hundreds other objects. 
Some had typical asteroidal and NEO orbits and reached Q<3 AU 
for several Myr. One object with initial orbit close to that of Comet 
88P/Howell after 40 Myr got Q<3.5 AU and moved in orbits with 
a«2.60-2.61 AU, 1.7< 9 <2.2 AU, 3.1<Q<3.5 AU, ewO.2-0.3, and i»5- 
10° for 650 Myr. If we consider this object, then for series n2 at d s <10 d 
the value of Td will be greater by a factor of 4 (i.e., «80 Kyr) than that 
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Table II. Mean probability P=10~ 6 P r of a collision of an object with a planet 
(Venus=V, Earth=E, Mars=M) during its lifetime, mean time T (in Kyr) during 
which q<a p i, T C =T/P (in Gyr), mean time Tj (in Kyr) spent in Jupiter-crossing 
orbits, mean time Td (in Kyr) spent in orbits with Q<4.2 AU, and ratio r of times 
spent in Apollo and Amor orbits. Results from BULSTO code at 10~ 9 <e<10 -8 
(marked as 1CT 9 ) and at £<KT 12 (marked as 1CT 12 ) and with RMVS3 code 
(Levison and Duncan, 1994) at integration step d s . In the case of asteroids, for 
the last four lines e o =0.05 and i =5°, and for other runs e o =0.15 and i o =10°. 
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in the corresponding line of Table II. The times spent by five specific 
objects that have large probabilities of collisions with the terrestrial 
planets while in IEO, Aten, A12 (l<a<2 AU, g<1.017 AU), Apollo, 
and Amor orbits are presented in Table III. 

With RMVS3 at d s <W days for 2P run, the value of P for Earth for 
one object presented in line 1 of Table III was greater by a factor of 30 
than for 250 other objects (see Table II). For series nl with RMVS3, 
the probability of a collision with Earth for one object with initial 
orbit close to that of Comet 44P/Reinmuth 2 was 88.3% of the total 
probability for 1200 objects from this series, and the total probability 
for 1198 objects was only 4%. This object (line 2 in Table III) was not 
included in Table II with 7V=1199 for nl. 

Table III. Times (in Myr) spent by five objects in various orbits and probabilities 
of their collisions with Venus (p v ), Earth (p e ), and Mars (p m ) during their lifetimes 
Th (in Myr). 
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For BULSTO at ee[10~ 9 -10~ 8 ] two objects (lines 3-4 in Table III) 
with the largest probabilities were not included in Table II for 2P at 
iV=501, and for 10P at iV=2149. The probabilities of collisions of these 
two objects with Earth and Venus (see Table III) were greater than for 
9350 other objects combined (0.17 for Earth and 0.15 for Venus). Large 
values of P for Mars in the nl runs with BULSTO were caused by a 
single object with a lifetime of 26 Myr. Ipatov (1995) obtained the 
migration of JCOs into IEO and Aten orbits using the approximate 
method of spheres of action for taking into account the gravitational 
interactions of bodies with planets. 

The times spent by 22000 JCOs in Earth-crossing orbits with a<2 
AU were due to a few tens of objects with high collision probabilities. 
With BULSTO at 10~ 9 <e<10~ 8 six and nine objects, respectively from 
10P and 2P series, moved into Apollo orbits with a<2 AU (A12 orbits) 
for at least 0.5 Myr each, and five of them remained in such orbits 
for more than 5 Myr each. The contribution of all the 9337 other 
objects to A12 orbits was smaller. Among the 9352 JCOs considered 
with BULSTO, only one and two JCOs reached IEO and Aten orbits, 
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respectively. Only one object in series n2 (line 5 in Table III) got A12 
orbits during more than 1 Myr. 

For the nl series of runs, while moving in JCO orbits, objects had 
orbital periods P a <20 yr (Jupiter-family comets) and 20<P a <200 yr 
(Halley-type comets) for 32% and 38% of Tj=0.12 Myr, respectively. 

Some former JCOs spent a long time in the 3:1 resonance with 
Jupiter and with 2<a<2.6 AU. Other objects reached Mars-crossing 
orbits for long times. So JCOs can supply bodies to the regions which 
are considered by many scientists (Bottke et al., 2002) to belong to the 
main sources of NEOs. The probabilities of collisions of bodies with 
the Earth per unit of time, i.e. the values of 1/T C , were of the same 
order for JCOs and resonant asteroids. Therefore, mean eccentricities 
and inclinations of Earth-crossers were similar for former TNOs and 
resonant asteroids. With BULSTO the mean probability of collisions 
with the Earth for the 5:2 resonance was 1/3 of that for the 3:1 reso- 
nance at e o =0.05 and this difference was greater by a factor of several 
at e o =0.15 (see Table II). 

The ratio Ps of the number of objects colliding with the Sun to the 
total number of escaped (collided or ejected) objects was less than 0.015 
for the simulations, except for Comet 2P/Encke, Comet 96P/Machholz 
1 from n2 series, and resonant asteroids. In the case of close encounters 
with the Sun, the values of Ps obtained by BULSTO and RMVS3 and 
at different e and d s were different, but all other results were similar, 
as probabilities of collisions of objects with the terrestrial planets were 
usually small after their close encounters with the Sun. 

The results presented in the paper were obtained for direct modelling 
of collisions with the Sun, but usually they are practically the same if 
we consider that objects disappear when perihelion distance q becomes 
less than the radius r$ of the Sun or even several such radii (i.e., we 
checked q<ksrs, where ks equals 0, 1, or another value). The only 
noticeable difference was for Comet 96P from n2 series and a smaller 
one for Comet 2P. For n2 series, several runs, in which there was an 
appreciable difference in time spent in orbits with Q<4.7 AU for ks=0 
and for ks=l (the times can differ by a factor of several), were not 
included in Tables 2 and 4. This difference was due to Comet 96P. 
Eccentricity and inclination of this comet are large, and they become 
even larger after close encounters with the Sun, so usually even for these 
runs the collision probabilities of objects with the terrestrial planets 
were not differed much (by more than 15%) at ks=0 and %=1. There 
were three runs, for each of which at ks=0 a body in orbit close to that 
of Comet 96P was responsible for 70-75% of collision probabilities with 
the Earth, and for ks=l a lifetime of such body was much less than for 
ks=0. Nevertheless, for all (~ 10 4 ) objects from n2 series, at ks=0 the 
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probabilities of collisions with the terrestrial planets were close to those 
at ks=l, even if we consider the above runs. The difference for times 
spent in Earth-crossing orbits is greater than that for the probabilities 
and is about 20%. For all runs at 2P series, the difference in time spent 
in orbits with Q<4.7 AU for ks=0 and for ks=l was less than 4%. In 
2P series of runs (and also for the 3:1 resonance with Jupiter), at ks=0 
we sometimes got orbits with z>90°, but practically there were no such 
orbits at ks>l (Ipatov and Mather, 2003a-b). For Comet 96P we found 
one object which also got i>90° for 3 Myr. Inclinations of other orbits 
initially close to the orbit of this comet did not exceed 90° . 



3. Trans-Neptunian Objects in Near-Earth Object Orbits 

Using the results of migration of TNOs obtained by Duncan et al. 
(1995), considering the total of 5 x 10 9 1-km TNOs with 30<a<50 AU 
(Jewit and Fernandez, 2001), and assuming that the mean time for 
a body to move in a Jupiter-crossing orbit is about 0.12 Myr, Ipatov 
(2001) found that about jV Jo =10 4 1-km former TNOs are now Jupiter- 
crossers, and 3000 are Jupiter-family comets. Using the total times 
spent by N simulated JCOs in various orbits, we obtained the following 
numbers of 1-km former TNOs now moving in several types of orbits: 

Table IV. Estimates of the number of 1-km former TNOs now moving in several 
types of orbits 
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3700 
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For example, the number of IEOs NiEOs=Nj tjEo/(Njtj), where 
tiEO is the total time during which Nj former JCOs moved in IEO 
orbits, and Njtj is the total time during which Nj JCOs moved in 
Jupiter-crossing orbits. The number of former TNOs in Apollo and 
Amor orbits can be estimated on the basis of nl and n2 runs. The 
number of NEOs with diameter d>l km is estimated to be about 1500 
(Rabinowitz et al., 1994) or 1000 (Morbidelli et al., 2002). Half of NEOs 
are Earth-crossers. Even if the number of Apollo objects is smaller by 
a factor of several than that based on nl and n2 runs, it is comparable 
to the real number (500-750) of 1-km Earth-crossing objects (half of 
them are in orbits with a<2 AU), although the latter number does not 
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include those in highly eccentric orbits. The portions of objects in Aten 
and A12 orbits are much greater in our 2P runs than in other runs. Our 
estimates of these portions are very approximate. The above estimates 
of the portion of former TNOs in NEO orbits are relatively large (up 
to tens of percents), but it is also possible that the number of TNOs 
migrating inside solar system could be smaller by a factor of several 
than it was earlier considered. 

Comets are estimated to be active for T aci ~10 3 -10 4 yr. T act is smaller 
for closer encounters with the Sun (Weissman et al., 2002), so for 
Comet 2P it is smaller than for other Jupiter- family comets. Some 
former comets can move for tens or even hundreds of Myr in NEO 
and asteroidal orbits, so the number of extinct comets can exceed the 
number of active comets by several orders of magnitude. The mean time 
spent by Encke-type objects in Earth-crossing orbits is >0.4 Myr. This 
time corresponds to >40-400 extinct comets of this type. Note that the 
diameter of Comet 2P is about 5-8 km (Fernandez et al., 2000; Lowry 
et al., 2003), so the number of 1-km Earth-crossing extinct comets can 
exceed 1000. The rate of a cometary object decoupling from the Jupiter 
vicinity and transferring to an NEO-like orbit can be increased by a 
factor of several due to nongravitational effects (Asher et al., 2001; 
Fernandez and Gallardo, 2002). 

Based on the collision probability P = 4x 10~ 6 we find that 1-km 
former TNOs collide with the Earth once in 3 Myr. This value of P 
is smaller than that for our nl and n2 runs and does not include the 
'champions' in collision probability. Using P = 4 x 10 -6 and assuming 
that the total mass of planetesimals that ever crossed Jupiter's orbit 
is ~lOOm0, where is the mass of the Earth (Ipatov, 1993), we 
concluded that the total mass of water delivered from the feeding zone 
of the giant planets to the Earth could be about the mass of Earth 
oceans. 

Our runs showed that if one observes former comets in NEO orbits, 
then most of them could have already moved in such orbits for millions 
of years. Some former comets that have moved in typical NEO orbits 
for millions or even hundreds of millions of years, and might have had 
multiple close encounters with the Sun, could have lost their mantles, 
which caused their low albedo, and so change their albedo (for most 
observed NEOs, the albedo is greater than that for comets; Fernandez 
et al., 2001) and would look like typical asteroids, or some of them 
could disintegrate into mini-comets and dust. 

Chen and Jewitt (1994) noted that while cometary splitting is some- 
times associated with the close passage of a comet by the Sun, it is 
also known to occur at heliocentric distances of up to 9 AU. At 10m, 
the near-Earth flux is more than two orders of magnitude greater than 
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power law extrapolated from larger sizes. Levison et al. (2002) obtained 
that majority of comets evolved inward from the Oort cloud must physi- 
cally disrupt, but Jupiter-family comets do not appear to disrupt at the 
same rate. Bailey (2002) consider that some long-period comets become 
innert and hence evolve into low-albedo objects resembling asteroids, 
and another alternative is that Oort cloud comets may easily break up 
into unobserved smaller bodies or dust. 

From measured albedos, Fernandez et al. (2001) concluded that the 
fraction of extinct comets among NEOs and unusual asteroids is signif- 
icant (>9%). Rickman et al. (2001) and Jewitt and Fernandez (2001) 
considered that dark spectral classes that might include the ex-comets 
are severely under-represented and comets played an important and 
perhaps even dominant role among all km-size Earth impactors. 

We conclude that the trans-Neptunian belt can provide a signif- 
icant portion of the Earth-crossing objects, or the number of TNOs 
migrating inside solar system could be smaller than it was earlier con- 
sidered, or most of 1-km former TNOs that had got NEO orbits dis- 
integrated into mini-comets and dust during a smaller part of their 
dynamical lifetimes if these lifetimes are not small. 

This work was supported by NRC (0158730), NASA (NAG5- 10776), 
INTAS (00-240), and RFBR (01-02-17540). We thank the anonymous 
referee for helpful remarks. 
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